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Metliod and device for measuring, determining and controlling 
flatness of a metal strip 

5 Technical Field of the Invention 

The invention relates to a method and a device for measiir- 
ing, determining and controlling flatness of a substant3_a.llY 
long and flat sheet or strip of material, such as copper-, 
10 steel, aluminium or other materials such as plastics, in a 
continuous or semi -continuous production process. 

More particularly it is a method and device for measuring, 
determining and controlling flatness in a cold or hot 
15 rolling mill- 

The invention further deals with a method and device for* 
measuring, determining and controlling flatness using hiaman- 
machine interfaces and control algorithms . 

20 

Background of the invention 

In the rolling of strip and sheet materials it is common 
25 practice to roll a material to the desired dimensions in a 
cold or hot rolling mill stand and then feed the resulti-ng 
strip to a coiler where it is wound up into a coil. In 
figure 1 a flatness control is illustrated that typicall-y 
comprises of three main parts, measurement, control and 
30 actuators, that correct the flatness error. The difficulty 
with the flatness control consists on achieving a very ^ood 
flatness in spite of certain known limitations. The reason 



wo 2005/064270 



PCT/SE2004/002059 



2 

why the flatness is not always perfect, is that all three 
parts have a number of limitations , 

The measurement is for example often discrete both across 
5 the width of the strip, as it has a limited number of 

measurement zones, and along the strip, as there are only a 
limited number of measurements done per rotation (say, four 
measuring zones in the rotational direction) . The control 
has a complex mathematics that cannot always follow at the 

10 speed of the strip and the actuators are only correcting 

certain type of errors and it is not sure they can correct 
any shape of flatness error. For example, a bending and a 
tilting actuator can correct only a shape that is a sum of a 
parable and a straight line, but of course the flatness 

15 error can in principle have any irregular shape. 

A metal strip that is subject to different degrees of 
reduction across its width will be elongated in varying 
length over different sections of the strip. Normally, the 
2 0 reduction results in certain pre-defined stress profiles for 
the strip. 

The difference between the measured strip flatness and a 
target flatness curve is defined as the ''flatness error". 
25 Practically, the flatness error of the strip can be gauged 
by a special measuring roll, e.g. a Stressometer roll, an 
ABB product- 
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Prior Art 

Several different methods have been developed to correct the 
flatness error by a so-called flatness control operation . 
5 All of these methods are based upon local modifications of 
the gap between rolls. 

The tension in the strip between a mill stand and a coiler 
is carefully monitored and it is known to measure tension 
10 distribution across a strip in order to regulate the 
flatness of the rolled material. 

In US 3,481,194 Sivilotti and Carlsson disclose a strip 
flatness sensor. It comprises a measuring roll over which 
15 the strip passes between a mill stand and, for this example, 
a coiler. The measuring roll detects the pressure in a 
strip at several points across the width of the strip. The 
pressure represents a measure of the tension in the strip. 
The measurements of tension in the strip result in a map of 

2 0 flatness in each of several zones 'across the width of the 

strip . 

US 4,400,957 discloses a strip or sheet mill in which 
tensile stress distribution is measured to characterize 
25 flatness. The measures of flatness are compared to a target 
flatness and a difference between measured flatness and 
target flatness is calculated, as a flatness error. The 
flatness error is fed back to a control unit of the mill 
stand, so as to regulate and control flatness in the strip 

3 0 in order to approach a zero flatness error. 
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Discrete flatness sensors and rolling mills with several 
actuators have drawbacks in that the sensors are not 
measuring continuously. This means that certain types of 
flatness errors are not detected. In addition to this, due 
5 to physical limitations, actuators cannot correct any type 
of flatness error, even if these are detected. 

One further disadvantage is increased complexity in sensor 
processing, actuator response prediction and the need of 
10 precise mathematical mappings between these sensors and 
actuators. This complexity can be overcome by an 
optimization operation of the unflatness and control. 

Normally the optimization of the flatness is done by 
15 different compensations of the space mapping function. These 
involve a great deal of numerical and process control 
methods. Typical examples are artificial neural networks, 
multivariable control, adaptive control, etc. These are 
difficult for plant operators to understand and handle. 

20 

Moreover, though additional information may be available for 
the operator from e.g. visual inspection of the strip, the 
operator cannot easily translate this information into the 
existing mathematical formalism of the flatness control. For 
25 example, the operator may see certain patterns on the strip 
like ''quarter-buckles" , but it is not obvious how he may use 
this information to amend the flatness control paradigm. 

In the following the background of the generic control 
3 0 algorithms for control equipments will be described for the 
purpose, of later understanding of the invention. In doing 
this, references will be made to the following parameters: 
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N = Number of measurement zones (number of measured samples 
along the width of the strip) 
P = Number of actuators 
5 Jbi = Flatness error for zone i (i=l,...,N) The vector b={jbi} is 
the flatness error space 

a.i^j = Actuator effect corresponding to measurement zone i 
for actuator j , , ..,,N; j=l , P) 

A = The mill matrix, with components ai,j used for space 
10 mapping 

jCj = actuator parameter to calculate, normalized in range 
[0, . . ,1] , the vector x={xi} is the actuator space 
A*^ = Means the transpose of a matrix A 

= It is the space of real n-dimension vectors {xi, X2,"-, 
15 Xn) . 

The term space denotes a vector space. We use the vector 
space of matrices of different sizes and the vector space of 
polynomials - 

20 

The mill matrix A is obtained from the mill manufacturer, 

from measurements or identification operations, and the 
flatness error matrix b is measured on-line with the 
flatness measurement device (such as a measurement roll) . 

2 5 These matrices are related to the actuator parameter matrix 

X by the following formula: 

Ax^b (1) 

3 0 Flatness control means generally finding the vector x when 

the vector A is a given constant and b is measured on-line. 
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Since A is not a square matrix (in general N > P) , there are 
more equations N than unknowns P for the system of equations 
represented by (1) . That means, in general, (a) there is no 
a unique solution and (b) the solution is not exact, i.e. 
5 there can be found an x such that b is "'close" to A:>c. Here 
""close" is in terms of a Euclidian distance expressed as a 
norm in the space - 

If jc varies and if the columns of A are linearly 
10 independent, then the product A?c describes a space on if 
called ""tiie column space of A'' or ""the range of A" - 
notation R(A) , or ""the span of the columns of A". The 
column space gives all the flatness error correction 
possibilities a rolling mill can have with its actuator 
15 equipment described by matrix A. It is known from elementary 
geometry in 2D/3D spaces and from functional analysis in N- 
space that the closest distance is given by an orthogonal 
projection: the smallest distance between b and Ax is the 
projection of b on the space of Ax, see figure 2. 

20 

Two arbitrary vectors v and w are orthogonal when their 
inner product denoted <v, w> is zero. For a matrix vector 
space, the inner product is defined as the product v'^w, so 
orthogonality means v'^w=0. For our case, let r = &-Axbe the 
2 5 residual (the error between a solution Ax and the vector b) . 
The smallest residual is the orthogonal projection of r onto 
Ax, which according to the definition is: 

0 = A^r = A'^(fe-^) (2) 

30 

By separating the terms, we obtain: 
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A^Ax^A^b (3) 

This is called in literature the noirmal equation. The factor 

a'^A eR^^ is called the Gram matrix of A (is always a 
5 symmetric square matrix) . Multiplying (2) with (a'^A)'^ we 
obtain x: 

x = <iA^Ar'A'b = A^''b (4) 

10 The term A^^ = (A^A)"^ is called the left pseudo-inverse 
matrix of A. 

The formula (4) above is used for instance for the "'model- 
based method" in the existing ABB Stressometer 5.0/6.0 

15 flatness control equipment and also by a majority of other 
producers of flatness control equipment. Other descriptions 
of this method appear under the names such as least square 
method", "'orthogonal vectors" or ""residual square 
minimization", however, all of these are less general 

2 0 formulations of the vector space and orthogonal projection 
described above. For example orthogonal vectors may be 
obtained by a Gram-Schmidt orthogonalization procedure; 
however the orthogonalization is just another way of 
computing the pseudo-inverse. The proposed method applies 

25 for all these computations or others that are equivalent. 

Image morphing (short for image metamorphosis) is a popular 
technique for creating a smooth transition between two 2D 
images. Morphing is the process of changing one object to 
30 another one. The method can be applied to volume-based 
representations of objects and the result is a smooth 
transition between two 3D objects (3D metamorphosis or 3D 
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morphing) . There are several different algorithms to create 
the transition depending on the intended function. 

If the objects have the same number of polygons and vertices 
5 and the morphing is done to 100 % the result is that the 
first object is transformed to an object identical to the 
one you morphed it with. I.e. if the two objects are a 

sphere and a cube you can morph the sphere to be identical 
to the cube if they have the same number of polygons and 
10 vertices. You can also choose to morph the objects to i.e. 
5 0 % and then you end up with a merge of a sphere and a 
cube . 

If the objects have different number of polygons or vertices 
15 you have to build a map between the vertices of one object 
and the vertices of the other. Therefore you do not end up 
with a morph shape looking identical to the one you morphed 
it with even if you morph the objects to 100 %. 

2 0 The method can be applied to volume-based representations of 

objects and the result is a smooth transition between two 3D 
objects (3D metamorphosis or 3D morphing) . There are several 
different algorithms to create the transition depending on 
the intended function. 

25 

Stunmary of tilie invention 

As strip products with a high degree of flatness are essen- 

3 0 tial for the competitiveness of mills today, the object of 

the present invention is to solve above indicated problems 
and present a new method and a device for improving flatness 
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control to achieve a very good flatness by measuring, 
determining and controlling the flatness of a substantially 
long and flat sheet or strip of material . 

5 Flatness Control, FC, essentially maps points moving in the 
multidimensional space of the measured flatness error into 
points determined in the multidimensional space of the 
actua.tors such that certain optimality conditions are met. 
This is the traditional mapping used today in FC, we call 
10 this as the ^ first mapping'. 

The first mapping is often difficult to optimize due to the 
nature of the actuator space and due to the sensor data vs. 
reference values used for the computation of the flatness 
15 error. For instance the complete action of all the actuators 
is not linear in the individual actuators, it cannot produce 
the optimal mapping or the mapping (possible combinations of 
actuator positions) is not unique. 

2 0 Another factor that limits the control is the discrete 

nature of the measurement. A discretization is inherent with 
most of the current measurement technologies both across the 
strip width (the measurement zones of the measuring roll) 
and along the strip length (sensors measuring when in 
25 contact with the strip) . 

The present invention proposes a solution to the problems 
discussed above by introducing a second mapping. This second 
mapping can be made from reference models such as video 

3 0 samples or mathematical functions to output space 

optimization matrices. The reference models enter into the 
compiatation of the main mapping. The result is increased 
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acc-uracy and locality of the flatness control. A further 
advantage is a vast increase in sampling granularity that 
translates into more accurate and timely flatness control. 

5 The method of the invention is achieved by, 

- that a second mapping is made to optimization algorithms 
acting in actuator model and sensor spaces using visual or 
image information, and 

- that the measured flatness vector is enhanced with 
10 information from a virtual strip. 

The present invention applies to any type of flatness 
measurement that is discrete in at least one dimension 
(along strip width, strip length or both) , such as 
15 measurement device or contactless measurement devices with 
discrete sampling. 

The present invention creates a mathematical map between the 
sensor space and the actuator space that improves 

2 0 substantially the deficiencies with known controlling 

methods. This comes from that extra information is entered 
into the system from the user, or from a visual sensor, 
about the nature of the flatness error and that this 
information is used in a better way to optimize the actuator 

25 behavior, for any type of known controller. Accordingly, a 
visual equipment and visual operation, on sensor data, such 
as sensor fusion or morphing is used in the method according 
to the invention, as well as a control architecture that 
implements a space mapping . 



30 
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A method accoording to the present invention will increase 
the control SE)eed substantially as the response time 
decreases . 

5 

Brief Descrlpt:;ion of Drawings 

For better uncier standing of the present invention, reference 
will be made to the below drawings / figures . 

10 

Figure 1 illustrates a known flatness control equipment - 
Figure 2 illustrates a simple orthogonal projection. 
Figure 3 illustrates a control system according to the 
invention. 

15 Figure 4 illustrates a mathematical model generated strip. 
Figure 5 illustrates a number of examples of flatness 
profiles. 

Figure 6 illustrates the fine-tuning of the reference model 
or the synchronization of the sampling. 

2 0 Figure 7 illustrates an example of a continuous strip 

prediction. 

Detailed Descoription of Preferred Embodiments 

25 

Flatness control, as used in metal rolling mill industry is, 
from a mathema-tical perspective, a mapping of an abstract 
space of sensor-originated vectors into a space of 
actuators. When fine-tuning the control to specific mill 

3 0 conditions, one uses traditional paradigms corresponding to 

the implemented mathematical level of the control, such as 
"transfer function", "poles and zeroes in the state space", 
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''dynamic behavior", ''least-squares" and different 
parameterization or penalty functions. With the advent of 
the IT and virtual reality architectures, it is increasingly 
the case tha.t virtual strip models are available from video 
5 captures or mathematical models . Thus the idea is to move 
the tuning operation from the control abstract level up to 
the real physical properties that can be seen by an operator 
when using different visualization and virtual reality 
techniques . 

10 

The present invention presents a method to use the 
information contained in the virtual strip models to improve 
the flatness control algorithm. The main features of the 
invention arre : 

15 

• The measiared flatness vector is enhanced with information 
from the virtual strip using a morphing technique 
(enhancement of the first space mapping) . 

• By this operation the sampling granularity increases. That 
20 result in. better and more accurate flatness control. 

• Since sampled data is available in shorter time interval, 
the contr-ol speed increases (response time decreases) 
substantially - 

• The virtual strip is mapped to a vector space optimization 
25 algorithm. This in effect is a second vector mapping that 

performs a fine, local adjustment of the control 
optimization (creation of a second space mapping) . 

In short, the invention states that increased flatness 
30 control precision can be achieved by a two-way space mapping 
in which visual or other continuous or near continuous a- 
priori strip information enhances (a) the measurement (b 



wo 2005/064270 



PCT/SE2004/002059 



13 



vector) and a traditional control loop and (b) the 
optimization algorithms that tells how x is computed such 
that to minimize the rezidual |Ax - b| . 

5 In figure 3 is illustrated a control system according to the 
invention, where: 

• A is a flatness measurement device, that produces a vector 
of flatness values per each measurement zone. 

• B is a Morphing lanit (software in existing equipment or 
10 dedicated device) , that is a part of the HMI . 

• C denotes reference strip models used for morphing. These 
are picture files generated from video or using laser 
captures or results from mathematical modeling. 

• D are target stress profiles (reference) . For the purpose 
15 of the method, the granularity of the target stress 

profile is the sa.me as the granularity of the reference 
strip models. 

• E is flatness control computer executing the algorithm 
described above. 

20 • F are space conversion matrices that describe 

mathematically the local properties of the strip (edge, 
middle etc) . Each matrix is mapped to a reference strip 
model . 

• G is a mill model, that is part of the existing flatness 
25 control, and it ±s a matrix (^'the Mill Matrix"). 

• H are actuators that correct flatness. 

• I is a device sensing the actual strip and converting this 
information into a format that allows comparison with 
reference strip models. This can be a measurement device 

3 0 such as visual camera, infrared camera or other devices 

capable of sensing strip relevant information, or a 
trained operator with visual access to the strip. 
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• J is a device that allows comparison between the actual 
strip sensed and converted according to I, and reference 
strip models. This operation, again, can be manual or 
automatic. However, in case of automatic operation, the 

5 comparison was done at the time of the design of the 

automatic device. 

• K is a device that selects the best match between the 
sensed actual strip and the reference strip models. 

10 The mill operator has extra information about the nature of 
the flatness fault, for example from a video camera or by 
visual inspection. By contrast, the measured value contains 
only a truncated part of all this information due to the 
discrete sampling. A video frame can contain information at 

15 sub-millimeter granularity while the flatness measurement 
may have only granularity between 10 cm and up to 0.25 m. 
For example a video camera (or visual inspection) can show 
clearly so-called ''quiarter buckles" but when sampling with a 
flatness measurement device, these buckles are not 

20 completely detectable. This is illustrated in figure 6 on a 
strip picture generated from a video capture, discrete roll 
measurement samples naay fall anywhere inside or outside 
buckles such that the buckles cannot be deduced from roll 
measurements alone. The visual information can be generated 

2 5 also from mathematical models (so called grrey, white or 
black models) , an example of a mathematically generated 
model strip is illustrated in figure 4. According to the 
invention all types of generated or photo /video captured 
images can be used. 

30 

According to the present invention, this information is 
added to the measured information from the reference model 
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in order to improve the space mapping optimization. The 
information is typically morphed. 

The operator has a number of reference models that have as 
5 output stress profiles, see figure 5, that can be morphed 
with the real measurement data. These models are obtained 
for example from digital translations of video captures. 

To accomplish the morphing the selected reference model need 
10 to be scaled and synchronized to match the measured data. 

One alternative to do this is to make a visual mark on the 
strip (from video capture) where the measurement was done. 
This can be done for example by .using a laser ray. 

15 

In the graphical user interface the operator (or a pattern 
matching algorithm) can fine-tune the reference model hori- 
zontally to match the measured data. This method is shown in 
figure 7 . 

20 

The morphing algorithm modifies the reference model 
according to the measurement data resulting in a data set 
that is closer to the reality. 

2 5 The morphing function can be calculated in different ways 
but essentially every measured data point affects several 
reference model points gradually. 

This new data set can be used to calculate a flatness error 
30 that is more accurate as compared to the existing system. 

The target profile is subtracted from the morphed profile (a 
vector subtraction) in order to compute a flatness error. 
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This operat±on can be performed if the granularity of the 
target and the morphed signal is the same. High granularity- 
target profiles can be obtained by simply filling in 
repetitively the same profile until the required granularity 
5 is achieved, or by using pictures or video captures from a 
de post- facto strip processing as described in the Patent EP 
1 110 63 5 Al Method and device for controlling- flatness 
(Jonsson, Meyer, 1999). 

The resulting flatness error has a much higher granularity 
10 than what is traditionally obtained from the original 

measured samples. This has the effect that the flatness 
control will be more accurate, faster and has a better 
prediction capacity. 

15 The mill operator has extra information about the nature of 
the flatness fault, e.g. from a video camera or visual 
inspection. By contrast, the measured value does not contain 
all this information due to sampling. For example a video 
camera (or visual inspection) can show so-called ^'quarter 

20 buckles" bixt when sampling with a discrete sampling 
measurement device, these buckles are not completely 
detectable . 

The operator has a number k of reference models of pre- 
25 defined strress profiles. These models are obtained from for 
example digital translations of video captures. 

Information about the nature of the flatness fault can be 
used to design the flatness control algorithm to concentrate 
30 on compensating zones where the problem is most significant. 
If for example the strip has the pattern ^^wave edges" the 
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actuators need to compensate strip edges rather than the 
center . 

The present invention proposes that there is established a 
5 one-to-one mapping: each reference model of the pre-defined 
stress profiles has associated a matrix that describes the 
weight of each measurement zone for the flatness error 
correction. These weight matrices are used in the flatness 
control algorithm. 

10 

The operator, or a pattern matching software, decides witch 
reference model is best to use. The reference model is 
mapped to one of the weight matrices. This mapping is the 
second vector space mapping-, and it is a key element of this 
1 5 invent i on . 

There are some patents from e.g. Hitachi / Siemens , that use 
certain pattern matching techniques. However, these patents 
are focused on the use of certain algorithms (fuzzy, 

2 0 neuronal nets) for enhancing measurement and are not 

referring to the control algorithm as this invention does. 
In other words, these patents claim improvements in what we 
call the first space mapping. By contrast, the present 
invention claims that substantially increased flatness 

2 5 control accuracy can be achieved by the dual space mapping. 

Mathematically, that means each reference model Mi (where i 
is the index of the model) , is mapped to a unique matrix K±. 
The matrix Ki is a square matrix of size N x N (N being the 
30 number of measurement zones) . The matrix has zero in all 
positions except the main diagonal, where there is the 
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weight desired for the action of actuators on the particular 
zone of size. 

The near- continuous space obtained after morphing, gives, 
5 after a measurement sampling, a degree of prediction since 
the continuous subspace between two samples is now 
available. The prediction horizon is ecjual to at least the 

size of the reference strip model. That means the control 
algorithm can be improved in this strip region, by 
10 performing an autoregressive control scheme. The operator 
has the visualization and command possibilities as shown 
schematically in figure 6 . 

The strip is generated on the operator station using 

15 computer graphics tools (shadowing, material texture, 

translations etc) , The strip is shown moving, in the same 
fashion as a video would show (so called virtual reality) . 
The input to this function is the measured data from the 
measuring roll and the morphed surface. The measuring 

20 sensors touch the strip e.g. in the positions denoted as A, 
B, C, D in figure 7. The reference model used for morphing 
has a length 1 . Figure 7 shows the case when the prediction 
horizon is 1. This length covers 3 measurement samples A, B 
and C. The colors within the prediction area are mapping the 

25 difference between the actual measurement and the predictive 
values, in the usual way done with loggers, Matlab programs 
etc. For example red color means large error, yellow means 
medium, green is low error and the colors have gradients in 
between. The mill operator, by modifying the prediction 

30 horizon in the range [0, I], and the morphed reference 

model, can verify which case is best suited for rolling. 
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To summarize, the invention involves the following; 

- a second mapping - from reference models to output space 
optimization matrices - virtual strip is directly mapped to 
a vector space optimization algorithm, 

5 - reference models enter into the computation of the main 
mapping - the measured flatness vector is enhanced with 
information from the virtual strip using a morphing 
technique, 

- use of reference pictures/models/descriptions of the strip 
10 to enhance flatness error - morphing technique, 

- use of reference pictures/models/descriptions of the strip 
with an associated weight matrix to concentrate on 
compensating zones where the problem is most significant, 

- space change . 

15 

The method according to the invention may, at least partly, 
be performed under control of a set of computer readable 
instructions or code means contained in a computer program 
storage device for making a computer or processor perform 
2 0 any of the steps of the above described method. 

The invention may also use a computer readable product for 
carrying out the method according to the invention , 

25 While the present invention has been described in terms of 
the preferred embodiments, the invention is not limited 
thereto, but can be embodied in various ways without 
departing from the principle of the invention as defined in 
the claims . 
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